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ABSTRACT 

Active galactic nuclei (AGN) drive fast winds in the interstellar medium of their 
host galaxies. It is commonly assumed that the high ambient densities and intense 
radiation fields in galactic nuclei imply short cooling times, thus making the outflows 
momentum-conserving. We show that cooling of high- velocity shocked winds in AGN 
is in fact inefficient in a wide range of circumstances, including conditions relevant to 
ultra-luminous infrared galaxies (ULIRGs), resulting in energy- conserving outflows. 
We further show that fast energy-conserving outflows can tolerate a large amount 
of mixing with cooler gas before radiative losses become important. For winds with 
initial velocity v m > 10, 000 km s _1 , as observed in ultra-violet and X-ray absorption, 
the shocked wind develops a two-temperature structure. While most of the thermal 
pressure support is provided by the protons, the cooling processes operate directly only 
on the electrons. This significantly slows down inverse Compton cooling, while free free 
cooling is negligible. Slower winds with u in ~ 1, 000 km s _1 , such as may be driven by 
radiation pressure on dust, can also experience energy-conserving phases but under 
more restrictive conditions. During the energy-conserving phase, the momentum flux 
of an outflow is boosted by a factor ~ v- m /2v s by work done by the hot post-shock gas, 
where v s is the velocity of the swept-up material. Energy-conserving outflows driven by 
fast AGN winds (y- m ~ 0.1c) may therefore explain the momentum fluxes P Lagn/c 
of galaxy-scale outflows recently measured in luminous quasars and ULIRGs. Shocked 
wind bubbles expanding normal to galactic disks may also explain the large-scale 
bipolar structures observed in some systems, including around the Galactic Center, 
and can produce significant radio, X-ray, and 7-ray emission. The analytic solutions 
presented here will inform implementations of AGN feedback in numerical simulations, 
which typically do not include all the important physics. 
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1 INTRODUCTION 

Recent observations show compelling evidence of galaxy- 
scale outflows driven by active galactic nuclei (AGN) . These 
include winds of ionized, neutral, and molecular gas in ultra- 



luminous infrared g 


alaxies (ULIRGs; 


Fischer et al. 


2010 


Feruglio et al.||2010 


Rupke & Veillcux 2011; Sturm et al. 


2011) and quasars ( 


Fu & Stockton||2009; 


Moe et al. 


2009 


Dunn et al. 2010; Bautista et al.|2010||Alexander et al. 


2010 


Villar-Martin et al.|2011 1. Possibly related fast outflows have 



2007 Tripp et al.|2011 \ and in Lyman break galaxies ( Hain 



line et al.|2011[ ). While quasar feedback has long been postu- 
lated to play important roles in galaxy formation, including 



establishing the Mbh — a relation (e.g., Silk & Rees 1998 
Wyithe fc Loeb||2003| |Di Matteo et al.||2005|> and trimcat- 



ing star formation in massive galaxies ( Springel et al.|2005 



Hopkins et al. 1120081), it is only recently that observations 



have directly revealed the impact of AGN outflows on galaxy 
scales. Furthermore, the data set on such outflows is poised 
to explode in the near future as the Herschel Space Observa- 
£onj]] continues its mission and the Atacama Large Millime- 
ter Arraij^ (ALMA) becomes operational. It is thus impor- 
tant to develop a theoretical understanding of the dynamics 
of AGN-driven winds. 

The most basic dynamical question regarding AGN 
winds is whether they conserve energy or momentum. Al- 
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though energy is conserved globally, the wind can in prin- 
ciple radiate away the thermal energy generated when it 
shocks with the surrounding interstellar medium (ISM). On 
the other hand, momentum cannot be radiated away. In 
the limit in which the shocked wind energy is rapidly ra- 
diated away, we term the outflow "momentum-conserving." 
If radiative losses are negligible, we say that the outflow is 
"energy-conserving." Real outflows may of course be some- 
what intermediate between these two limits. Nevertheless, 
the degree to which AGN outflows conserve energy has very 
significant implications. In particular, the momentum flux 
of the material swept up in energy-conserving outflows in- 
creases with time owing to work done by hot shocked gas. 
An analogous phenomenon operates in supernova remnants 
(SNRs), in which the momentum of the remnant at the end 
of the Sedov- Taylor (energy-conserving) phase can exceed 
the momentum of the explosion ejecta by a factor as high as 
Cioffi et al.|1988 l. Such momentum boosts have 



50 (e.| 



a large effect on the efficiency of stellar feedback in galaxies 
(e.g., |Ostriker fc Shetty|2011| |Hopkins et al.|20l"2| ). 

Both observations and theoretical modeling suggest 
that galaxy-scale AGN outflows have momentum fluxes well 
in excess of the momentum flux output radiatively by the 
central black hole (BH), P ra d = Lagn/c. Observationally, 
the "momentum boost" P/P la d ranges from ~ 2 to ~ 30 



in local ULIRGs dominated by AGN (Rupke & Veilleux 
2011 Sturm et al.|20lT l and in luminous quasars for which 
the outflow properties have been measured accurately ( Moe 



et al.|2009||Dunn et al.|2010[[Bautista et al.|20T0||Faucher 



Giguere et al.|20TT I. The existing measurements are summa- 
rized in greater detail in §3.4| Although the uncertainties are 
large, collectively these measurements indicate that AGN- 
driven, galaxy-scale outflows may commonly have momen- 
tum fluxes ~ 10Lagn/c. On the theoretical side, simulations 
that model AGN feedback by depositing radial momentum 
on scales ~ 100 pc have found that comparably high mo- 
mentum fluxes are needed to reproduce the normalization 
of the M BH - a relation ( |DeBuhr et~aT1|2011b|a"| ). If AGN 
outflows are energy-conserving, it is possible that they are 
radiatively launched with momentum flux ^ Lagn /c in the 
nucleus, then boosted to ~ 10Lagn/c on galaxy scales in a 
Sedov- Taylor-like phase. 



The original work of Silk & Rees ( 1998 1 on quasar 



feedback implicitly assumed that the outflow was energy- 
conserving. Subsequently, King ( 2003 1 argued that in- 
verse Compton cooling by the quasar radiation field 
should efficiently cool the shocked wind, making the out- 
flow momentum-conserving out to ~ 1 kpc, but energy- 
conserving past this radius (|King et al. 2011 Zubovas 



fc King||20l2| ). |Silk fc Nusser| ( |2010[ ) argued that energy- 
conserving outflows are not possible in galaxy bulges, al- 
though the cooling times in their numerical models appear 
too short in comparison with the analytic theory developed 
herein. In this work, we focus primarily on fast (> 10, 000 
km s _1 ) nuclear winds and show that cooling is less impor- 
tant than previously appreciated. 

Such winds may be launched from accretion disks 
around central black holes and give rise to ultra-violet 
(UV) broad absorption lines (BALs; Murray et al.|[T995[ ). 



Observed BALs have typical maximum velocity offsets 
10, 000-30, 000 km s _1 , and even higher velocity absorption 
has been detected in some cases (e.g., Weymann et al.|1981 



Turnshek|1988l|Gibson~t al. 2009). However, we will refer to 
fast nuclear winds generically in this paper since the winds 
need not be directly connected to an accretion disk (for ex- 
ample, if driven by radiation pressure on a larger scale). 
Furthermore, fast winds need not manifest themselves in 
the form of classical BALs. X-ray observations also find ev- 
idence for AGN winds with velocities ~ 0.1c, albeit with 
lower significance and with more ambiguous interpretation 



(e.g.,|Chartas et al. 12002] |Pounds et al.|2003]|Tombesi et al 



20101. We also briefly consider slower winds (~ 1,000 km 
s~ ), as may be accelerated by radiation pressure on dust 
(e-g. 



Roth et al.|2012 \ . These winds too can experience sig- 



nificant energy-conserving phases, but under more restric- 
tive conditions. 

The cooling of the shocked wind in AGN involves im- 
portant conceptual issues which have not been previously 
discussed in this context. Furthermore, key physical pro- 
cesses (in particular, inverse Compton cooling and two- 
temperature plasma effects) are not included in many ex- 
isting numerical simulations of AGN feedback. Simulations 
of isolated galaxies or galaxy mergers, and most acutely cos- 
mological simulations, also typically do not resolve all the 
relevant scales, including the wind shock. Thus, we focus 
here on exploring the basic physics of galaxy-scale outflows 
analytically. This work will inform the implementation of 
AGN winds in future numerical simulations, which will be 
best suited to address potential three-dimensional effects in 
more realistic settings. 

The plan of this paper is as follows. In Sj2] we investigate 
the cooling physics of fast shocked winds, and show that the 
weak collisional coupling between protons and electrons sup- 
presses the inverse Compton cooling rate significantly. We 
also show how mixing of the shocked wind with cooler gas 
modifies the effective cooling rate. In Sj3] we use the cooling 
rates derived in fj2] to calculate representative outflow solu- 
tions. We show that the AGN-driven outflows are expected 
to be energy-conserving for a wide range of parameters rele- 
vant for observed systems. We argue that this conclusion is 
robust to the effects of mixing, and show how only modest 
confinement of the hot gas is necessary to produce momen- 
tum boosts comparable to those measured. We conclude in 
SQwith a discussion of the implications of our analysis, in 
particular with respect to recent observations of galaxy-scale 
AGN winds. 

A series of appendices support the main text. Appendix 
|A"| presents a general derivation of the conditions for energy 
and momentum conservation for self-similar outflow solu- 
tions. Appendix[B]summarizes predictions for emission from 
shocked wind bubbles. Finally, Appendix [C] outlines how to 
incorporate the physics of AGN winds in numerical simula- 
tions. 



2 COOLING OF THE SHOCKED WIND 

Our first goal is to derive the cooling rate of a fast, shocked 
AGN wind. 



2.1 Problem setup and shock structure 

We focus first on the case of fast nuclear winds with launch- 
ing speeds Vi n > 10,000 km s" 1 . For concreteness, we ana- 
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than all others in the problem. For instance, the accretion 



shocked 



QSO^ 




shocked 
wind 



sw 



Figure 1. Schematic illustration of the outflow structure. The 
AGN wind is launched in the galactic nucleus with velocity u; n . 
The wind is shocked at R BW (velocity v BVJ ). A second shock, at 
radius R s (velocity v s ), is driven into the ambient ISM. The two 
shocks are separated by a contact discontinuity at R c . The swept- 
up ambient material piles up in a shell at R s . However, it is the 
cooling of the shocked wind (not the shocked ambient medium) 
that determines whether the outflow conserves energy or momen- 
tum. 



lyze the following idealized problem. At t — 0, a central black 
hole with radiative luminosity Lagn turns on and drives a 
constant, spherically-symmetric wind with initial velocity 
Vin and mass outflow rate Mi n . For reference, the Eddington 
luminosity (LEdd) is related to the black hole mass (Mbh) 
via 

47tGMbHTO p C _.-,„.. -,„46 -1 ( M BH 



1.3 x 10 46 erg s" 1 



10 s M Q 



(1) 

where G is Newton's constant, m p is the mass of the proton, 
c is the speed of light, and or is the Thomson cross sec- 
tion. We also define the Eddington accretion rate MEdd = 
£Edd/(^c 2 ), where 77 = 0.1 is the radiative efficiency. The 
momentum flux of the small-scale wind is parameterized by 

Lagn (2) 



If Lagn = ^Edd, then 

Tin = V 



(3) 



M Bd d 

is related to the mass loading of the wind. 

The ambient, pre-shock medium is also assumed to 
be spherically symmetric, with a radial gas density profile 
p s (R). The details of the small-scale wind acceleration are 
unspecified, but could owe to radiation pressure on spectral 
lines, free electrons, or dust, as well as hydro-magnetic forces 
(we return to the acceleration mechanisms in [4.31. We as- 
sume that this acceleration occurs on scales much smaller 



disk wind model of Murray et al. ( 1995 1 predicts that the 
wind is launched at a radius ~ 0.01 pc. Since the inter- 
esting physics in our calculations occurs at radii where the 
escape velocity is <C v- m , we can also neglect gravity in our 
first-order calculations; the numerical calculations of [|3]will 
include a more accurate treatment. 

Before proceeding, it is important to understand the 
structure of the outflow that develops. Studies of the analo- 
gous problem of fast stellar winds revealed that the outflow 
is characterized by both a forward and a reverse shock, sep- 
arated by a contact discontinuity (see Fig. [I] e.g., Weaver 



et al. 19771. The outermost shock is the forward shock, 



which propagates in the ambient medium. The innermost, 
reverse shock is the shock first encountered by the gas ac- 
celerated by the BH; it is also termed the wind shock. We 
denote the radius of the wind shock by R svl , the radius of the 
contact discontinuity by _R C , and the radius of the ambient 
medium shock by R s . Jl sw and 7? a have velocities v sw and v s , 
respectively, with respect to the BH. As the outflow prop- 
agates outward, it sweeps up ambient ISM. This results in 
a shell of swept-up material at m R a and causes the bubble 
to slow down its expansion. Except at very early times, we 
thus have v B <C v m and iw <C UinnThen, the velocity of the 
wind shock relative to the wind itself is v ws = Vi n — v sv , w Wj n . 

Neglecting for the moment two-temperature effects, 
which we demonstrate are important in the next section, 
the temperature of the post-shock gas depends sensitively 
on the shock velocity. For a strong shock of velocity v sn in 
a monatomic gas, the Rankine-Hugoniot jump conditions 
imply 



Tsh(«sh) = Yok mpV!ih 



1.2 x 10 10 K 



t'sh 



30, 000 km s" 



(4) 

where fi is the mean molecular weight. In general, v s <C 
Vi n and the shocked wind is much hotter than the shocked 
ambient medium. For instance, while the swept-up material 
may decelerate to a few 100 km s _1 (or less), the velocity 
of the wind shock remains nearly constant at v ws ~ Win > 
10, 000 km s~ . This is key to the question of whether AGN 
outflows conserve energy or momentum, because it is much 
harder to cool the T > 10 9 K shocked wind than the cooler 
shocked ambient medium. In particular, it is the cooling of 
the shocked wind, rather than the shocked ambient medium, 
that determines whether the outflow is energy or momentum- 
conserving (see also King et al.||20li |. Whether or not the 
shocked ambient medium cools has no significant impact on 
the global outflow dynamics, and for our purposes R c w _R S 
(e.g., Koo & McKcc 1992a I. We thus make the identification 
i? c — > Re in the following. 

To determine the cooling rate of the shocked wind, we 
need its characteristic density, temperature, and radius from 
the BH. Although the exact structure of the wind bubble 
is non-trivial, the analysis is considerably simplified, yet re- 
tains the essential physics, by adopting the following approx- 
imation. Since the sound crossing time of the shocked wind 



3 The last inequality may be transiently violated if the shocked 
wind quickly loses its pressure support and R Bvi — > R B on a time 
scale as short as ~ (R s — R S w)/vi n , but this does not affect our 
arguments. 
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(t s ~ [Rs — Ravr]/vin) is short, the shocked wind pressure is 
nearly uniform between f? sw and R B - Furthermore, |Koo &:| 
McKee ( 1992b I show that the gas density p is approximately 
constant in the intermediate region Rs-w R "^C j a re- 



sult supported by the numerical calculations of Falle ( 1975 1 



Briefly, this follows because mass conservation and adia- 
baticity imply that R 2 v is constant in that region. Since the 
flow is steady, pR 2 v is also constant, implying that p must 
be constant. As the pressure P oc pT is uniform, T is approx- 
imately uniform as well. The properties of the shocked wind 
are thus adequately described by constants p sw (i), T aw (t), 
and P sw (i) , where the density and temperature are given by 
the shock jump conditions at 7? S w In particular, 



= 4p„(i? sw ) = 



M in 



(5) 



Before cooling of the shocked wind becomes important, the 
thermal energy of the wind bubble is conserved, E b = 
(l/2)Li n t, where we assume that 1/2 of the energy goes 
into bulk motion of the swept up shell, 1/2 goes into ther- 
mal pressure of the shocked wind bubbler] and L- m = 
(l/2)M in vf n . The thermal pressure of the shocked wind is re- 
lated to E h via P sw = 2E h /3V h , where V h = 47r(-R?-i?| w )/3. 
Balancing the ram pressure from the unshocked wind with 
the thermal pressure of the shocked wind bubble at i? sw , we 
obtain 



Rsw — 



Lin 

E b v h 



(Rs — 7?s 



1/2 



(6) 



When _R SW <C Rs, which applies when the wind shock is not 
radiative, _R SW ~ R s (v s /v in ) 1/2 (using t ~ Re/v e ). 

Finally, the cooling of a Lagrangian element of shocked 
wind occurs principally at ~ Rs, so this is the relevant radius 
to evaluate the inverse Compton cooling time below. 



2.2 Cooling of the two-temperature plasma 

The key issue in determining the outflow type is whether the 
shocked wind cools. We distinguish between three regimes: 
energy-conserving, momentum-conserving, and the interme- 
diate partially radiative bubble stage. The limit realized de- 
pends on the ordering of three time scales: the flow time 
t n ow = Rs/vs, the crossing time t cr = R B v,/v ln , and the cool- 
ing time of the shocked wind t coo i (Koo & McKee 1992b I. 
If icooi 3> £ n0 w > tor, then the outflow is energy conserv- 
ing. If tcooi <C tcr < tflow, the outflow conserves momen- 



tum. In the intermediate partially radiative bubble stage 
tcr < tcooi < tflow, most of the shocked wind cools and is 
compressed in a thin shell at Rs. However the wind shock 
itself is non-radiative and most of the volume is filled by the 
portion of the wind that has not had time to cool yet. In this 
stage, the outflow conserves neither momentum nor energy. 

For fast winds, free free emission and inverse Compton 
scattering in the AGN radiation field are the most important 
cooling processes (e.g., King 2003; Ciotti & Ostriker 1997, 
2001 also emphasized the importance of Compton scattering 



4 This is a simple approximation to more exact calculations; 
|Weaver et al.| ( |1977^ show that for adiabatic bubbles in a uni- 



for AGN inflows and outflows). In practice, free free cool- 
ing is subdominant at the high temperatures T > 10 9 K of 
shocked winds with Vi n > 10, 000 km s~\ but we include it 
in our calculations for completeness. 

A complication neglected so far in the context of AGN 
winds concerns the coupling between protons and electrons 
in the shocked wind. This is important because most of the 
cooling occurs via the electrons, but the protons carry most 
of the thermal pressure. Observations of supernova remnants 
in fact indicate that the electrons and protons can be signif- 



icantly decoupled after fast shocks (e.g., Ghavamian et al. 
20071, as is also the case in the solar wind (e.g., Schwartz 



et al. 19881. If the protons in shocked AGN winds are not 



well coupled to the electrons, then their thermal energy can 
be trapped in the wind, in which case the outflow effectively 
conserves energy. Such a situation is analogous to low ra- 
diative efficiency accretion flows, in which two-temperature 
plasmas are invoked as a means to prevent the energy from 



being radiated away (e.g., Rees et al. 1982 Narayan & Yi 
1995||Quata"ert||2001| ) 



We first discuss in more detail the empirical and the- 
oretical evidence regarding two-temperature effects. Some 
measurements of Balmer-dominated SNR shocks indicate 
that T e can be as low as the minimal heating case ~ 
(m e /m p )T p for fast shocks with Vsh 2,000 km s _1 (e.g., 
Ghavamian ct al. 2007). However, the evidence is ambigu- 
ous for SNR shocks in general. For example, observations of 
the reverse shock in SN 1987A suggest an electron-proton 
temperature ratio as high as T c /T p « 0.14 — 0.35 ( |France| 
et al. 2011 1. In the solar wind, the faster shocks are measured 



to heat electrons to less than 10% of the proton tempera- 
ture (Schwartz et al. 19881. Markevitch| ( |2006| ) presented ev- 
idence for electron heating faster than predicted by Coulomb 
collisions in the bow shock of the "Bullet" galaxy clus- 
ter. However, the Mach number M w 3 of the bow shock 
is much lower than expected for AGN winds. While elec- 
tron heating is not well understood theoretically, particle- 
in-cell simulation demonstrate that electrons can be quickly 
heated to T e > 0.1T P by plasma instabilities in collisionless 
shocks (Riquelme & Spitkovsky 2011 1. Importantly, though, 
the same simulations show that magnetic turbulence decays 
rapidly in the post-shock region (IChang et al. 20081. As 



form medium E b K, (5/ll)Li n t. 



we show below, our results do not depend significantly on 
the level of "prompt" electron heating at the shock, but 
rather on the assumption that Coulomb collisions dominate 
afterward. We also note that even if some residual collision- 
less heating persisted well past the shock itself, the ratios 
T p /T e ~ 10 that we find below assuming Coulomb collisions 
only is comparable to (or less than) the ratios observation- 
ally determined in most fast shocks in SNRs and in the solar 
wind. Our calculations are thus relatively robust to uncer- 
tainties in collisionless processes and are well motivated by 
the available theoretical and empirical data. 

Since most of the pre-shock kinetic energy is carried 
by the protons, minimal electron heating implies that those 
particles initially receive the majority of the shock heat, 
T p « 2T a h(t'in), with electrons initially heated to T c ~ 
Tp(m e /m p ). As mentioned above, it is possible that plasma 
instabilities near the shock heat the electrons to higher 
temperatures. Our arguments are however not sensitive to 
the electron temperature immediately past the shock, as 
long as coupling in the post-shock region is dominated by 
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Coulomb collisions. This is shown explicitly below, in Fig- 
ure]^ where we explore the case of complete heating at the 
shock, T p = T c . Assuming that Coulomb collisions in fact 
dominate in the post-shock region, electrons and protons 
achieve equipartition (T c ~ T p ) on a time scale 



tei 



3m m p 
8(2 7 r) 1 /2 npe 4 1nA 



fcXl kT^ 



3/2 



where 



In A w 39 + In 



(i0 10 k) 2 ln (lcm- 3 ) 



(7) 



(8) 



( Spitzer|1962"| >. 

Equation |7| assumes a fully ionized plasma consist- 
ing of protons and electrons only, an approximation which 
we adopt throughout. We fiducially assume In A = 40 in 
this work. In our treatment of the shock, we further assume 
that there is a magnetic field sufficiently strong to suppress 
thermal conduction, but small enough to have a negligible 
dynamical impact. This allows us to treat electron-ion equi- 
librium locally. In this section, we use expressions appropri- 
ate in the limit of non-relativistic particles, which allow us 
to derive simple analytic approximations. It is not a priori 
clear that this is a valid assumption since, for example, the 
electrons could be heated to relativistic speeds in reality. We 
have implemented an alternative set of equations using the 
expressions from Narayan & Yi (1995), which are valid in 



both the relativistic and non-relativistic regimes, and con- 
firmed that the results obtained here are accurate. 
The thermal free free emissivity is 

6 ff (T c ) = 1.4 x lQ~ 27 T c 1/2 Cnlg B , (9) 

where all terms are in cgs units, g~B is the average Gaunt 
factor, and C is a clumping factor (which we both set to 1 
fiducially). The rate of energy transfer between the electrons 
and photons owing to Compton scattering is 



dllph ( kT e 

at \m c c z 



(10) 



5 kT c 2 fcT x 

1 + o 2 _ i71 " 2 

2 m e c z m c c^ 



where T x Compton temperature of the radiation field (e.g., 
Sazonov & Sunyaev 2001) and 

Lack 



(7 ph = 



(11) 



4tyR 2 c 

is the radiation energy density. 

Let U e and Up be the electron and proton energy den- 
sities, respectively. For simplicity, we assume that the elec- 
tron and proton distributions are individually thermal at all 
times, so that U a = 3n c kT B /2 and Up = 3n p kT p /2 (neutral- 
ity implies that n p = n c ). Assuming that the plasma is of 
uniform and constant density, the temperature evolution of 
the plasma is described by the following system of equations: 



dT a 
dt 



3kn c 



dU p h 



dt 



+ es(T e 



+ 



dT p 
dt 



T c — T p 

tei 



(12) 



(13) 



To develop an understanding of how the shocked wind 
cooling proceeds, we consider an idealized model in which 
T p = 2T sh (v in ), T c = Tp(m c /m p ) at t = 0. 



Figure [2] shows the evolution of the proton and electron 
temperatures for the representative cases of Lagn = 10 46 
erg s _1 (a Mbh = 10 8 Mq black hole radiating near the 
Eddington limit), v in = 10, 000, 30, 000 and 50,000 km s" 1 , 
and Tin = 1. The different sets of curves in the top panel 
correspond to radii R B — 1, 10, 100, and 1,000 pc, from left 
to right. A Compton temperature T x = 2 x 10 7 K is assumed 
based on the average AGN spectrum from Sa zonov et al.| 
( |2004[ > . For these calculations, the density of the shocked 
wind is approximated by 



M in 



M in 



irR^vinmp nKivsirip 



Tin L AGN 

tyRsVs Vin crrip 



(14) 



The second equality assumes that the outflow is energy- 
conserving (i? S w ~ (vs/vin) 1 / 2 R s ). The calculations also as- 
sume an outer shock velocity v s = 1, 000 km s _1 . This is not 
necessarily self-consistent, as a given ambient density profile 
implies a relation between R s and v s ; self-consistent models 
are constructed in !j3] 

For reference, the vertical dotted lines show the rela- 
tivistic inverse Compton cooling time evaluated assuming 
equipartition (T c = T Bn (vin)) in each case, 



R _ 3m e c 

^Comp — 



32<T T [/phfcTe 

2 x 10 7 yr 



(15) 



Rs \ 2 / Lagn 



1 kpc / VlO 46 er S s_1 



30, 000 km : 



This is the relevant cooling limit neglecting two-temper ature 
effects for u in > 10,000 km s" 1 (e.g., as in |King|2003 |. 

Particularly for v in = 30,000 km s" 1 and v in = 50,000 
km s _1 , the protons take significantly longer than icomp 1° 
lose half of their thermal energy. Moreover, the effective cool- 
ing time of the protons increases with increasing Vi n , a trend 
opposite to standard inverse Compton cooling for relativistic 



electrons (eq. ( 15 1). The electrons in fact remain slow enough 



that their inverse Compton scattering time scale follows the 
non-relativistic expression to a good approximation. In con- 
trast, assuming T c ~ T p would imply a root mean square 
(rms) electron velocity {v 2 } 1/2 w A3(v 2 } 1/2 (where {v p } 1/2 
is the proton rms velocity), putting the electrons into the 
relativistic limit for these shock velocities. Thus, the effec- 
tive cooling time of the shocked protons is not the standard 
inverse Compton cooling time, but is longer owing to two- 
temperature plasma effects. 

We now derive an analytic approximation to the effec- 
tive cooling time of the shocked protons. First, it is appar- 
ent from Figure [2] that the electron temperature reaches a 
plateau before the protons cool significantly. This tempo- 
rary "equilibrium" electron temperature, Tc q , is set by a 
balance between Coulomb heating by the protons and in- 
verse Compton cooling in the AGN radiation field. Setting 
dT c /dt — 0, neglecting free free cooling, approximating in- 
verse Compton cooling by dll p h/dt = (4:kT e /m e c)aTn B Uph 
(exact for a thermal, non-relativistic electron distribution 
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Figure 2. Proton (solid) and electron (dashed) temperatures be- 
hind the wind shock, as a function of time since passing through 
the shock. Fiducial AGN parameters (Lagn = 10 46 erg s _1 , 
Tin = l)i a swept-up shell velocity v B = 1,000 km s , and 
a Compton temperature Tx = 2 X 10 7 K are assumed in all 
cases. The different colors in the top panel correspond to dif- 
ferent radii of the swept-up shell (R B , assumed fixed), and the 
three panels show the results for different wind shock velocities 
v V s ~ fin- The shocked wind number density is given by equation 
j!4| . For reference, the vertical dotted lines show the relativis- 
tic inverse Compton cooling time evaluated assuming equiparti- 
tion (T e = T sn (ni n )) in each case. The vertical dash-dotted lines 
show the effective proton cooling time taking into account two- 
temperature plasma effects (eq. |l7j). This effective proton cool- 
ing time is significantly longer than the standard inverse Compton 
cooling time, especially at high v ln when the electrons would be 
relativistic in the absence of two-temperature effects. Most curves 
assume minimal electron heating at the shock (T c = (m c / m p )T p ) . 
The thinner curves in the bottom panel show an example assum- 
ing full electron heating at the shock (T p = T e ), but Coulomb 
coupling only afterward. The effective proton cooling time is not 
sensitive to the initial electron temperature, as long as Coulomb 
collisions dominate past the shock. 



and T x « T c ), and in the limit T c cq /m c > T p /m p , 

1/5 



V(mA)XT p 2 



(16) 



,2/5 



1/5 



3 282/1 A\22 2 

m^fi e c (In A) Vi n r in 



2„,2 2 "I 1/5 



a^k 5 m p v B 



2.9 x 10 9 K ( ^ 
30 



2/5 
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^40~ 



2/5 
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The intermediate equality holds assuming U p h and n p are 
given by equations (111 (with R = R s ) and (14 1, respec- 



tively, and T p = 2T sh (v in ). 

This analytic approximation is in excellent agreement 
with the numerical results in Figure [2] (for the v in = 10, 000 
km s _1 case, proton cooling is significant before the plateau 
is reached) . Coulomb collisions with the electrons is the only 
process through which protons lose energy. The proton cool- 
ing time is therefore given by t PiCOO i w t c i(T° q ), i.e., 

4/5 r m >^ 3 D s 2 <C 8 ' lVS 



p,cool 



3^ / 7T\ 

8 V2/ 



1.4 x 10 a yr 



a|e 8 

R s 
1 kpc 

2/5 



(hlA)Vn 



Lagn 



(17) 
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Although this time scale derives from two-temperature ef- 
fects, it is nearly independent of the electron temperature 
immediately past the shock. This is because T° q is set by 
a balance between inverse Compton cooling and Coulomb 
heating of the electrons in the post-shock region. In partic- 
ular, independent of the initial temperature of the electrons 
at the shock, Tg q approaches the value given by equation 
( 16 1 in the downstream region (see Fig. [5J bottom panel) 
and thus the estimate of the effective proton cooling time in 



equation (17 1 applies 



Finally, consider synchrotron cooling. For a given elec- 
tron population, the ratio of the synchrotron to inverse 
Compton power is 

-Psynch _ Ub 



(18) 



^Comp 

where U B = B 2 /S TV is the magnetic energy density. For an 
AGN radiation energy density following equation (111, 



-fsynch 
-fComp 



B 2 R 2 c 

2I/AGN 

1.4 x 10" 



(19) 



B 

1 



R 



1 kpc y 

/ Lagn 
V 10 46 erg s- 



In ordinary spiral galaxies, B ~ 1 — 10 fiG (Beck 2001 1 and 
synchrotron cooling is clearly subdominant. In starbursts, 
the magnetic field can reach B ~ 1 mG (Thompson et al. 
|2006[ ). Compact starbursts, such as local ULIRGs, however 
have gas disk scale heights ~ 50 pc (e.g., Downes & Solomon 



19981. Provided that the magnetic field scale height is less 



than a few times this value, synchrotron cooling should also 
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be less efficient than inverse Compton cooling in those ob- 
jects. We thus neglect synchrotron cooling of the shocked 
wind. 



2.3 Mixing and other cooling limits 



Equation ( 17 1 is the cooling time appropriate in the absence 



of mixing with ambient gas. Here, we study the effects of 
mixing of cold gas with the hot shocked wind in the limit in 
which the mixed cold and hot gas is rapidly homogenized. 
Then, the amount of mixing can be parameterized by the 
fraction 



In 



M sw + Mcoid 
M RW 



(20) 



where M co id is the total mass of cold gas that is mixed with 
the shocked wind, and M sw = Mint. M co id can consist of 
contributions from cold clumps ablated by the wind and 
from mixing of the swept-up ambient medium owing to hy- 
drodynamical instabilities. 

Consider first a wind with Vi n > 10,000 km s" 1 which 
initially develops a two-temperature structure past the wind 
shock. Initially, the pressure support P sw tx n p T p is pre- 
served when mixing occurs. As shown in §2.2| T,f q also de- 
pends only on the product n p T p , and is thus also unaffected 
by mixing. In particular, in the absence of efficient colli- 
sionless coupling in the shocked wind, the electrons remain 
non-relativistic. We term this the Coulomb limit. 

On the other hand, the effective proton cooling time in 
the two-temperature plasma approximation, 

tp.cooi w ioi(T c cq ) oc — oc (21) 

rip /mix 

is reduced linearly with the amount of mixing. When the 
time scale for energy exchange between the electrons and the 
protons becomes shorter than the inverse Compton scatter- 
ing time, the latter process becomes the rate limiting step. 
We term this the inverse Compton cooling limit. It is realized 
when the non-relativistic inverse Compton cooling time 

37r m c c 2 _Ra 
2 <tt£agn 

R s \ 2 / £agn 

.ripe/ v 1046 er s s_ 

corresponding to 



NR 
Comp 



(22) 
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47r l/5 2 4/5 
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(23) 
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Since the tcSnp is independent of density and temperature, 
it is constant as mixing proceeds. 

As the amount of mixing becomes larger, / m i x — > oo, 
the increase in the shocked wind bubble density and the 
decrease in its interior temperature both shorten the free 
free cooling time, 

1/2 



3n c kT c /2 



,7 x 10° yr 
BbC 



To 
10 6 K 



n e 
1 cm -3 



(24) 



Multiplying equation ( 14 1 by / m i x to estimate the shocked 



wind density after mixing (assuming energy conservation) 

and using T c = T s h(i>in)//mix, we solve for / m i x such that 
+ _ ,nr . 

Iff — 'Comp- 



J mix 



26.4 



(Cg B Tin) 2 / 3 \ 1,000 km s 



2/3 



(25) 



30,000 km s" 1 



4/3 



-3/2 
/ mix 



In this free free limit, ts oc / ~ 

It is important to note that the / m i x thresholds derived 
above determine the rate limiting step for cooling but do not 
imply that catastrophic cooling occurs. We show in [J3]that 
substantially more mixing is needed for the shocked wind to 
radiate away its energy in a flow time in realistic conditions. 

For slower winds, /^j x ap and / mix can be < 1. This 
simply implies that the cooling of the shocked bubble may 
never be in the Coulomb or inverse Compton limits. For 
example, if both I^^ p < 1 and / mix < 1, then the shocked 
wind cooling is dominated by free free emission even in the 
absence of mixing. This is the case for v m ~ 1, 000 km s _1 . 
For Vin < 1, 000 km s _1 , bound- free emission dominates over 
free free; we do not consider this limit further here. 

At each radius, there is a physical upper limit to /mix, 
realized when all the swept-up gas mass has been mixed 
with shocked wind. Let M s be the swept-up gas mass at 
-Rs, defined as the gas mass originally enclosed within that 
radius. Then, the maximum mixing parameter at R s is 



final/ r> 
J mix K^s 



Msw 



(26) 



In our numerical integrations (q3|, the mixing parameter at 
each radius is thus taken to be the minimum of the pre- 
scribed value and /mbf (-R s ). 



3 OUTFLOW SOLUTIONS 

We now use the results of the previous section on the cool- 
ing rate of the shocked wind to calculate outflow solutions. 
We focus here on numerical examples with parameters rep- 
resentative of AGN systems, but a general analysis of self- 
similar solutions in Appendix [^demonstrates that our con- 
clusions regarding energy conservation apply more generi- 
cally. In what follows, the gas density profile of the ambient 
medium is parameterized by 



Ps(R) = Po 



(27) 



We use Ro = 100 pc for all the numerical calculations in 
this work. 



3.1 Thin shell approximation 



Our approach is similar to that of Weaver et al. ( 19771, but 



we include additional cooling processes, gravity, and approx- 
imate the effects of mixing following §2.3| We assume spher- 
ical symmetry and that conduction is negligible. 
The basic equation of motion is 



d 
di 



(MsWs) =47Ti? s 2 (Pb-Po) 



GM s M t 

Ri 1 



(28) 



© 0000 RAS, MNRAS 000, 000-000 



8 Faucher-Giguere & Quataert 



where Pb is the thermal pressure in the shocked wind bubble, 
Po is the pressure of the ambient medium, M s = J^" dVp s is 
the small of the swept-up gas shell, and M t = MBH + M ga i(< 
R s ) accounts for the total gravitational mass within 7? a , in- 
cluding the central massive black hole and the surrounding 
galaxy. The gravitational potential of the galaxy is assumed 
to follow an isothermal sphere with velocity dispersion a, 



M gal (< R) 



2a l R 
G ' 



(29) 



In the calculations that follow, we fix this total potential 
even when varying the properties of the ambient gas density 
profile (eq. ( |27[ )). We do this because, as we argue in the next 
section, the effective gas density swept up by the outflow is 
not necessarily directly tied to the total gravitational mass. 
The velocity dispersion is set to a = 200 km s _1 , appropriate 
for Mbh = 10 s Mq black holes on the Mbh — <J relation 
(e.g., Giiltekin et al. 20091. This allows us to explore the 



dependences on the properties of the ambient medium while 
keeping the gravitational potential fixed. We have verified 
that instead assuming a total potential dominated by the 
gas mass does not change our results significantly. 

The thermal energy in the shocked wind region, E b , is 
related to its pressure, Pb, by 



E b = 2TrP b (Ri~Rl 



(30) 



The location of the wind shock, R sv , is given by pressure 
balance as in equation (|6j) . The rate of change of the energy 
in the shocked wind is determined by the energy injection 
rate, work done on the swept-up ISM, and cooling losses: 



E b = -M in < - 47vRtP h R s - L b , 



(31) 



where L b is the radiative cooling rate. The appropriate cool- 
ing limit depends on the amount of mixing, parameterized 
by / m i x ( £ 2.3). We define it as Z/2T, L-Comp, or Lb, depend- 
ing on whether / mix < f°™», < / mix < /«, 

fn 



> /mix, where: 



/iE b 



p,cool 



J Comp 



I^E h 

j-NR > 
Tomp 



L& 



^E b 
tff 



(32) 



At each radius, / mix is limited by /j£bf(-Rs) (eq. p6|). We 
do not keep track of the electron and proton temperatures 
separately, but instead use the approximation in equation 
( 17 I for £ Pl cool- We will also explore cases in which we neglect 
two-temperature effects on the cooling rate of t he s hocked 
wind. For those, i p , C ooi is replaced by tc om p ( ec L- |l5| ). Note 
that the mean molecular weight pre-factors are included in 
equation ( |32[ ) because the cooling time scales were defined 
relative to the thermal energy carried by the electrons (or 
protons) only. 



3.2 Escape along paths of least resistance 

For a galaxy-scale outflow with momentum P — 
Tl/AGNifiow/c (the normalization r is to be distinguished 
from the parameter n n characterizing the nuclear wind), 
momentum conservation implies that the average enclosed 
mass density is, assuming spherical symmetry and neglect- 



expansion along path 
of least resistance 
t"h~8 cm" 3 
(N H ~10 22 cnr 2 ) 



fi s =1 kpc 
^5=1,000 km/s 
P~15/. AGN /c 




swept-up ISM 
broad CO (-10%) 



Figure 3. Example of an outflow escaping preferentially along 
paths of least resistance above and below the dense molecular 
disk in a ULIRG system. Even though the gas density along the 
directions of escape is < 10 — 3 that of the molecular disk, the 
momentum boost can be boosted by factor ~ 15 during the ex- 
pansion to a kiloparscc. 



ing gravity, 



3r 
4ttR? cm, 



L-agn 

7 , - 



(33) 



0.5 cm 3 r 



L 



AGN 



Rs 



10 46 erg s- 1 J \1 kpc 

v, 

x 



1,000 km s- 1 



Here, we use nu to denote the total hydrogen number den- 
sity, even if it is in molecular form. In ULIRGs with ob- 
served outflows, this is much lower than the mean densi- 
ties nH ~ 10 3 — 10 4 cm" 3 of the dense molecu lar disks on 
scales ~ 500 pc (e.g., |Downes fc Solomon|l998^ p|Thus, the 
observed galaxy-scale outflows in these systems must have 
propagated along paths of least resistance. In the prototypi- 
cal example Mrk 231, the neutral outflow extends to Rs ~ 3 
kpc and has v s ~ 1,000 km s" 1 (iRupke & Veilleux||2011 1. 



As this system is viewed nearly face-on, the outflow is most 
likely observed as it escapes approximately normal to the 
dense molecular disk. 

To illustrate how little confinement is necessary to pro- 
duce a significant momentum boost, consider the example 
shown in Figure [3] based on the fiducial AGN parameters 
used previously. In this case, the shocked wind escapes pre- 
dominantly above and below a dense molecular disk. The 
molecular disks in ULIRGs have scale radii Rd ~ 500 pc 
and scale heights /id ~ 50 pc (e.g., Downes & Solomon 1998 1. 



The inward gravitational force on such a disk is 



d.grav 



GM d M t (< R d ) 



(34) 



Assuming that the AGN outflow is isotropic before encoun- 
tering the disk, the outward force it exerts is 



Pd,, 



(35) 



5 The dense molecular disks can extend to a radii > 1 kpc. The 
host galaxies (including the stellar disks and tidal features) are 



typically significantly more extended. 
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Figure 4. Dependence of the outflow solution on njj o (density at _Rq 



parameters, motivated by local ULIRGs, are: Lagn = 10 erg s , f n 



100 pc), a (where njj oc K 1 
: = 1, a = 1, ran,o = 100 cm" 



n , and / m i x . The fiducial 
in = 30,000 km s -1 , and 

T; n = 1- I n au cases, the total gravitational potential is described by an isothermal sphere with a = 200 km s — l , but the effective ambient 
gas density profile is varied. The solid line segments indicate periods during which the outflow is energy-conserving (t coo i > !fj ow ), the 
dashed segments indicate periods during which the outflow is in a partially radiative bubble stage (t cr < t coo i < £ now ), and the dotted 
segments indicate periods during which the outflow is momentum-conserving, fn each case, the curve starts at the free expansion radius, 
Rf (see AppendixjAj. The thin curves in the top left panel illustrate the impact of neglecting two-temperature effects on inverse Compton 
cooling. 



© 0000 RAS, MNRAS 000, 000-000 



10 Faucher- Giguere & Quataert 



where h^/Rd is the fraction of the momentum flux inter- 
cepted by the disk. Thus, 



~ T 



hdRd 



c GM d M t (<R d ) 
Lagn 



(36) 



0.07r 



10 46 erg s- 1 J V50 pcy V500 pc 



( Rd 



10 9 M 



Mt(< Rg) 
5 x 10 9 M Q 



Therefore, AGN outflows generally cannot eject the massive 
molecular disks in those systems in an instantaneous sense, 
unless r 3> 1. This is consistent with CO observations of 
ULIRGs, in which the outflowing gas with v s ~ 1,000 km 
s" 1 contributes only a small fraction ~ 10% of the integrated 
CO luminosity (Fcruglio et al.|2010 |, most of the flux being 
concentrated in a narrow component identified with the disk. 
On the other hand, r ^> 1 can be achieved if there is suffi- 



cient confinement of the hot gas. Indeed, Greene et al. ( 2011 1 



observe significant non-gravitational disturbances over the 
entire host galaxies of luminous obscured quasars at z < 0.5. 
Regardless, the wind is generally expected to preferentially 
escape normal to the disk. 

Even though the mean ambient density normal to the 
disk is < 10 -3 of the molecular disk density in this example, 
the kiloparsec-scale outfl ow still reaches a momentum flux 
P ~ 15-Lagn/c (see 3.4 below). Note that the implied col- 
umn density Nu ~ R B nn ~ 10 22 cm -3 is relatively modest, 
especially if we consider the column before feedback pro- 
vided by the wind exposes the AGN (e.g., Risaliti et al 
1999 Treister et al.|[2~009 i. In reality, the outflow propaga- 
tion depends on the details of the multiphase structure of 
the ambient ISM, and this aspect will be best addressed us- 
ing three-dimensional simulations. Nevertheless, this demon- 
strates that escape along paths of least resistance does not 
preclude large momentum boosts, unless these paths are ex- 
tremely under-dense. 



3.3 Representative examples 

In Figure [4] we illustrate the dependence of thin shell out- 
flow solutions on wh .o (density at Ro = 100 pc), a (where 
7iH oc R~ a ; eq. (27l), t>i n , and / m i x . The fiducial parame- 



ters for these calculations, motivated by local ULIRGs, are: 
£agn = !0 46 erg s" 1 , / mix = 1, a = 1, n H ,o = 100 cm -3 , 
Win = 30,000 km s _1 , and r; n = 1. The figure, in particu- 
lar, shows how maintaining a shell velocity v s « 1,000 km 
s" 1 at _R S « 1 kpc requires a relatively tenuous ambient 
medium, 71h(100 pc) « 10 cm -3 . In that case, the momen- 
tum boost (quantified by M E v s / (fm-^AGN /c)) is ~ 10, in 
good agreement with measurements of outflows in AGN- 
dominated ULIRGs (|Feruglio et al.|2010||Rupke fc Veilleux 



20TT)|Sturm et al.|2011 | 

In the top left panel (varying nn,o), we include exam- 
ples of the same calculations, but using the inverse Compton 



cooling time for single-temperature gas (eq. (15l) instead of 
the approximation for a two-temperature plasma (eq. ( IT ) ) . 



Without accounting for the slow down of inverse Compton 
cooling owing to two-temperature effects, the effective cool- 
ing time of shocked wind is underestimated, leading to solu- 
tions that depart from energy conservation for high ambient 



gas densities. This effect is potentially very important for 
outflows that are launched in the nuclei of ULIRGs. 

The examples in the bottom right panel show that mix- 
ing is not effective in removing the thermal pressure support 
of the shocked wind, with a discernible effect on the out- 
flow dynamics only for / m i x > 10 3 , and then only at radii 
R s > 1 kpc. Although mixing accelerates the cooling rate of 
the shocked wind via more rapid Coulomb collisions, two- 
temperature effects nevertheless play a role in keeping the 
electrons non-relativistic and slowing down inverse Comp- 



ton cooling ({2.3 1. In these examples, the large amount of 



mixing necessary to have a dynamical impact owes princi- 
pally to the fact that even when the cooling rate limiting 
step is free free emission, the pre-mixing temperature of the 
shocked wind is so high (T B h(vi u ) ~ 10 10 K) that mixing 
must reduce it by a large factor before tg/ta < 1. The im- 
pact of mixing is also limited at small radii by the constraint 
that the effective / m i x cannot exceed the value /.™t x (^s) at- 
tained when all the swept-up gas mass is mixed with the 



shocked wind (£ 2.3 1 . 

The lower left panel of Figure [4] includes examples of 
winds with Vi n = 1,000 km s _1 and Vi n — 3,000 km s _1 . 
These slower winds are proxies for outflows accelerated by 



radiation pressure on dust (e.g., Roth et al. 20121. While 



these examples initially experience momentum-conserving 
and partially radiative bubble phases, they eventually tran- 
sition to the energy-conserving regime and their momentum 
flux is also boosted significantly. On the other hand, these 
slower winds quickly slow down to velocities <C 1,000 km 
s _1 as they sweep up the ambient medium, a point to which 
we return in §4.3| 

Although it is difficult for an AGN outflow to instanta- 



neously unbind the bulk of the molecular disk (eq. (36 1), it 
remains an open question whether such outflows will even- 
tually empty the disk of most of its gas, as suggested by 
the short gas depletion time scales inferred by |Sturm et al.| 
( |2011[ ). Two limits are envisioned. In the first limit, the out- 
flow recently ejected the gas normal to the molecular disk, in 
which it was originally confined. In an instantaneous sense, 
the outflow has a large momentum flux and is highly mass- 
loaded. However, the wind subsequently escapes easily along 
paths of least resistance now cleared and has little global im- 
pact on the galaxy. In the other limit, the AGN wind is still 
only effective at blowing out the ambient gas in its imme- 
diate vicinity, but will eventually incorporate most of the 
ISM of the galaxy as gas is driven into the nucleus by gravi- 
tational torques. Such gravitational torques are expected in 
local ULIRGs, which were likely triggered by recent galaxy 
mergers (e.g., |Sanders et al.|1988"l|Mihos fc Hernquist|1994l ). 

We plan to investigate the global impact of AGN-driven 
outflows in more detail in future work. 



3.4 Momentum boost 

For outflows that are energy-conserving during most of their 
expansion, the expected momentum boost can be estimated 
simply using 



M s vi 



(37) 



again assuming that w 1/2 of the small-scale wind energy 
goes into bulk motion of the swept-up gas. Defining the 
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Figure 5. Momentum flux of the galaxy-scale outflow, relative 
to Lagn/c, as a function of velocity of peak absorption for sys- 
tems in which the outflow is believed to be driven by the AGN. 
Solid squares: Molecular outflows observed in OH by | Sturm et] 
[ah] (2011) , for the ULIRGs with strong AGN contribution. An 
uncertainty of a factor 4 is assumed in each case, dominated 
by the uncertainty on the mass outflow rate (the momentum 
flux is estimated as the mass outflow rate times the velocity of 
peak absorption). From left to right: IRAS 13120-5453, Mrk 231, 
and IRAS 08572+3915. Solid triangles: FeLoBAL quasars inter- 
preted with the radiative shock mod el of |Faucher-Giguere et ah] 
| |2011|, using data f rom|Moe et al.| ( |2009} , |Dunn et al.| l |2010| , 
and |Bautista et al.| ( |2010[ l. We assume a factor of 2 uncertainty 
from the outflow covering factor and the photoionization model- 
ing. From left to right: QSO 2359-1241, SDSS J0318-0600, and 
SDSS J0838+2955. The black solid curve shows the prediction 
P/(Lagn/c) = "in/21'peak implied by e nerg y conservation, as- 
suming V[ n = 0.1c and Pi n = Lagn/c (eq.l38J|. 



large-scale momentum flux as P s = M B v B and P m = ALnfin, 



Pin 



1 Wis 

2 v e 



(38) 



(see also Zubovas & King 2012| . Figure [5] compares this 
prediction with the observed outflows, as a function of the 
inferred velocity of peak absorption « poa k. The theoretical 
curve assumes Win = 0.1c. Interestingly, this prediction of 
energy conservation is in good agreement with the data 
(Veilleux & Rupke 2011, Sturm et al. 2011), both in normal- 
ization and in slope. As the numerical examples in Figure [4] 
show, gravity acts to suppress the momentum flux relative 



to the simple estimate in equation (381 when it becomes 
important as the flow decelerates. 



4 DISCUSSION 
4.1 Summary 

We have studied the conditions under which galaxy-scale 
AGN outflows conserve energy versus momentum. We fo- 
cused primarily on winds with initial velocity v- ln > 10, 000 
km s~ , as may be launched from the black hole accre- 
tion disk, but also considered slower examples as proxies 



for winds accelerated by radiation pressure on dust. Our 
study was motivated in part by observations of galaxy-scale 
outflows driven by AGN suggesting that they carry momen- 
tum fluxes many times the value P ~ Lagn/c expected if 
the outflows are accelerated by radiation pressure and each 
photon from the AGN scatters once with the surrounding 
medium (e.g., Sturm et al.|20lT l. Numerical simulations of 
AGN feedback also appear to require P ^> Lagn /c in order 
to reproduce the observed A/bh — o relation ( DeBuhr et al. 
2011b|a| >. If an outflow conserves energy rather than mo- 



mentum, then hot shocked gas can do work on the swept-up 
material, thus boosting the momentum flux, increasing the 
dynamical importance of the AGN outflow, and potentially 
providing a simple explanation for the observations of high 
values of P/ (Lagn/c). 

We first developed a theory of fast AGN wind bubbles 
under the assumption of spherical symmetry, building on 
analogous work in the context of stellar wind bubbles (e.g., 

~ 1992b[ ). The principal 



Weaver et al. 1977 Koo & McKee 



improvement of our analysis over previous work (e.g., King 
200H|King et al.|2011 l is a more complete treatment of the 
cooling of the shocked wind (§2p . We showed that the time 



scale for Coulomb coupling between protons and electrons 
in the shocked wind can be sufficiently long that the post- 
shock gas develops a two-temperature plasma structure, sig- 
nificantly modifying its effective cooling rate. The competi- 
tion between inverse Compton cooling and Coulomb heating 
of the electrons causes them to remain non-relativistic, and 
introduces Coulomb coupling as the rate limiting step for 
removing the thermal pressure support provided by the pro- 
tons. Both of these effects slow down the cooling of the pro- 
tons in the shocked wind relative to simpler estimates that 
assume that the electrons and the protons have the same 
temperature (and thus that the electrons are relativistic). 
We also carried out our analysis for ambient density profiles 
of arbitrary slope, rather than focusing on the isothermal 
case as in most previous studies. 

From our spherically-symmetric analysis, we conclude 
that fast AGN outflows conserve energy in a significantly 
broader range of circumstances than previously appreciated 
(ij3|. For instance, a wind with v ln = 30,000 km s" 1 and 
Pin = Lagn/c (n n — 1) driven by a black hole of mass 
Mbh = 10 8 Mq radiating at its Eddington limit is predicted 
to be effectively energy-conserving everywhere for ambient 
media of uniform density as high as uh ~ 8 x 10 6 cm -3 
(Appendix [Aj. If the outflow is energy-conserving, then the 
momentum boost of the swept-up material relative to the 
nuclear wind can reach P s /Pi n ~ Vi n /2v s . Thus, an energy- 
conserving outflow of velocity v B ~ 1,000 km s _1 driven 
by a nuclear wind with Vi n ~ 30,000 km s" 1 would have 
experienced a momentum flux boost by a factor up to ~ 15. 
Interestingly, these estimates are in good agreement with 
observations of galaxy-scale outflows believed to be driven 
by AGN in local U LIRGs (see Fig. [5] [Rupke fc Veilleux]20lT 
Sturm et aLp)TTj ) 



We also considered how three-dimensional effects might 
affect our conclusion. Mixing of cold gas with the shocked 
wind, either from ablation of clumps in the ambient medium 
or via instabilities at the interface with the swept-up shell, 
can accelerate the cooling rate. Depending on the amount 
of mixing, the rate limiting step for the shocked wind cool- 
ing can be either: 1) Coulomb collisions; 2) non-relativistic 
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inverse Compton cooling; or 3) free free cooling. Although 
non-relativistic inverse Compton cooling can be somewhat 
more rapid than the Coulomb-limited regime for outflows 
with Wi n 3> 10,000 km s _1 , it is still relatively slow; catas- 
trophic cooling can only occur when the cooling becomes 
dominated by free free emission. However, a large amount 



of mixing can be tolerated before this happens (J 2.3 1. For 
the parameters of our fiducial numerical calculations in i|3j 
motived by local ULIRG observations, a cold gas mass > lCr 
times that of the shocked wind mass must be mixed so that 
cooling is accelerated enough to have a discernible dynamical 
impact. This is much more than for the analogous problem 



of stellar wind bubbles (e.g., McKee et al. 19841 owing to 
the higher temperature of the shocked wind for the AGN 
case with Vi n ~ 0.1c, which requires more mixing before the 
cooling time is reduced to less than t cr or t now . 

Propagation of the outflow along paths of least resis- 
tance in an anisotropic, multiphase ambient medium ( "leak- 
age") tends to limit the momentum boost that can be 
achieved. However, only a modest amount of confinement is 
necessary to provide momentum boosts comparable to those 
inferred in ULIRGs and luminous quasars. In fact, spatially- 
resolved observations of the luminous quasar Mrk 231 show 
evidence of an AGN-driven outflow wit h v s ~ 1,000 km 
s _1 extending ~ 3 kpc from the nucleus (Rupke & Veilleux 



20111. The mean density along the path of propagation of 



this outflow must be much less than the density of the mas- 
sive molecular disk in the same system, since otherwise the 
outflow would have decelerated to much lower velocities (see 
eq. ( |33| l). As Mrk 231 is viewed nearly face-on, we are likely 
seeing the outflow as it leaks out normal to the galactic 
plane. Nevertheless, if the outflow is powered by a fast nu- 
clear wind, it has decelerated sufficiently to have experienced 
a momentum boost comparable to the inferred one. 

4.2 Importance of physical effects and comparison 
with previous work 

In this work we have studied several key physical effects 
previously neglected in models of the interaction between 
AGN winds and gas in their host galaxy. We now summarize 
each of these important effects in turn. 

The two-temperature cooling effects sl ow d own inverse 
Compton cooling for Vi n > 10, 000 km s" 1 ({ 2.2 1. Relative to 



the assumption that electrons and protons share the same 
temperature at all times, the cooling time of the shocked 
wind can be longer by a factor > 10 — 100 (Fig. [2J . Absent 
substantial mixing of the shocked wind with cooler gas, in- 
verse Compton scattering off the AGN radiation field (with 
proton cooling limited by Coulomb collisions with electrons) 
is the dominant cooling process. The two-temperature ef- 
fects thus directly affect the conditions for energy versus 
momentum conservation. In particular, nuclear winds with 
Win ~ 0.1c can be in the energy-conserving limit even for 
the high ambient gas densities and intense radiation fields 
characteristic of the centers of ULIRGs. The longer effective 
cooling time of the protons owing to two-temperature effects 
is robust to uncertainties in collisionless processes that may 
heat the electrons faster than Coulomb collisions near the 
shock. 

Previous analytic studies have focused primarily on 
isothermal ambient gas density profiles (a — 2). We have 



considered more general profiles with a < 2 and have shown 
that the isothermal profile is in fact a singular case. This 
is because it does not have a finite free expansion radius 
R{, so that outflow solutions retain memory of the initial 
conditions on large scales. With an isothermal sphere den- 
sity profile and a constant gas mass fraction / g , momentum 
conservation implies a constant swept-up shell velocity 
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,200 km s- 1 , 
such solutions are not con- 
> WmQ The smal l veloc- 



(King 2003). B ecause Ri — } oo, 
sistent with the assumption Vi n 3> v n 

ity v s — Vm implies a longer t now = R E /v B . Thus, indepen- 
dent of two-temperature cooling effects, we find that cool- 
ing is less important (than, e.g., King 2003) when the non 
self-similarity of realistic models is included. The numerical 
models in Figure [4] in fact show that energy-conserving so- 
lutions are obtained for parameters that reproduce observed 
galaxy-scale outflows even when neglecting two-temperature 
effects. 

In a real galaxy, the cold gas is flattened into a disk and 
the ISM is highly inhomogeneous. Thus, there is no clear 
reason to favor a = 2 as representative of the paths of least 
resistance along which galactic winds effectively propagate, 
even if the rotation curve is flat (j|3.2[|. For example, |Berger| 



et al. (20121 infer the pre-outflow gas density profile around 



the tidal disruption event Swift J164449. 3+573451 (assumed 
to originate from the central black hole) to have a slope 
a « 1.5 within 0.5 pc, with significant flattening farther 
out. Our general self-similar analysis in Appendix [A] shows 
that the transition from momentum to energy conservation, 
and vice versa, is sensitive to a. 

Mixing of the shocked wind with cooler gas tends to 
increase its cooling rate. However, the numerical examples 
in Figure [4] show that our conclusions are robust to a large 
amount of mixing (see also [ 2.3 1. 

The M B h - a relation (|Ferrarese fc Merritt|20"00"{ |Geb- 
|hardt et al.||2000) |Tremaine et al.||2002[ ) provides a con- 
straint on whether AGN winds conserve energy or momen- 
tum. King ( 2003 1 showed that the observed relation can 



be reproduced in a simple spherically-symmetric model of 
momentum-conserving outflowsQln an analogous model as- 
suming energy conservation, the predicted black hole mass 
lies below the Mbh — cr relation by more than one order 
of magnitude. Because in real galaxies the gas inflows that 
fuel central black holes are likely to occur through disks, 
isotropic black hole winds transfer only a fraction of their 
outward momentum flux to the inflowing gas. Thus, winds 



6 Self-consistent physical solutions do exist, but they differ from 
the self-similar expression in equation ( |39[ t . In particular, they are 
sensitive to the radius at which the wind is launched and start 
w ith v B = v ln 

7 |Murray et al.|p005| also showed that the observed Mbh-t can 
be derived from a simple Eddington limit argument for radiation 
pressure on dust. That argument however does not rely on the 
existence of a wind within the dust sublimation radius. 
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that are purely momentum-conserving may actually be in- 
sufficient to regulate black hole growth at the level implied 
by the Mbh — a relation. The estimate in equation ( 36 1 in- 
dicates that a momentum boost r ~ 10 — 15, such as can be 
attained in an energy-conserving outflow, may be necessary 
to stop inflows from the molecular disks in ULIRGs to the 
central black holes. Comparable momentum flux enhance- 
ments are in fact needed to reproduce the Mbh — a relation 
in three-dimensional simulations (e.g., DeBuhr et al.|2011b |. 

Finally, we note that although our analysis implies that 
winds driven by AGN are energy-conserving in a substan- 
tially wider range of conditions than found by King and 
collaborators, those authors also predicted a transition to 
energy conservation on scales > 1 kpc relevant for galactic 



winds (King et al. 2011 1. In addition, significant momentum- 
conserving phases are not ruled out by our arguments. If 
strong observational evidence for momentum conservation is 
found, our study of the conditions under which it is realized 
would provide useful constraints on the physical conditions 
in galactic nuclei. 

4.3 Driving mechanisms 

Our models do not depend on how the wind is launched 
on small scales, but we did focus on the case of high initial 



velocities, v- m > 10, 000 km 



Furthermore, in order to 



explain galaxy-scale momentum fluxes P > 10Lagn/c in- 
ferred from observations by hot gas confinement, we require 
that the wind starts with P[ n ~ Lagn/c. We now discuss 
the evidence that such winds are realized and the physical 
mechanisms that can drive them. 



UV BALs are observed in up to 40% of quasars (e.g., Dai 
et al. 2008 ) and could be present in the majority of quasars 
if this fraction represents a viewing angle effect ( jWeymann 
et al. 19911. The physical properties of the absorbing gas 



are generally not well known. Photoionization modeling of a 
subset of low-ionization absorbers indicates that the kinetic 
luminosity of the outflowing gas can reach a few percent of 



Lagn (Moe et al. 2009 Dunn et al. 2010 Bautista et al 



|2010[ ). While the absorption in those cases likely arises from 
the interaction of the quasar outflow with the ISM of the 
galaxy ( Faucher-Giguere et al.|2011 1 , rather than in the im- 
mediate vicinity of the black hole, the energetics do suggest 
that AGN can drive dynamically important winds. The elec- 
tron scattering optical depth around black holes accreting 
near the Eddington limit is furthermore naturally r cs ~ 1 
(e.g., King & Pounds 2003). Since the rate of gas inflow from 
galaxy to nuclear scales is not tuned to the Eddington limit, 
T e a may in fact exceed unity by a significant factor in many 
cases. X-ray observations confirm that this is the case in the 



prototypical example Mrk 231 (Braito et al. 2004). If the nu- 
clear wind is accelerated by radiation pressure on electrons, 
its momentum flux on small scales is thus P[ n > Lagn/c in 
these circumstances. The assumptions of our analysis may 
thus be generically realized in luminous quasars. 

Dust in the central regions of ULIRGs is also gener- 
ally optically thick to the UV photons that dominate the 
radiative output of AGN, as well as to reprocessed infrared 
radiation (e.g., IScovillc 2003 Thompson et aL"||2005[). Like 



the electron scattering opacity, dust opacity is continuous in 
wavelength and can thus scatter the bulk of the AGN lu- 
minosity. Furthermore, the large optical depths in ULIRGs 



imply that the momentum flux imparted on the surround- 
ing medium can potentially be several times Lagn / c owing 
to photon trapping (Roth et al. 2012; but see Krumholz & 
Thompson 2012, Novak et al. 2012). Thus, realistic simula- 
tions of AGN feedback in galaxies should include radiation 
pressure on dust in addition to the fast nuclear winds ana- 
lyzed in this work. 

On the other hand, radiation pressure on dust cannot 
readily explain AGN winds commonly observed in UV and 
X-ray absorption with velocities 2> 1,000 km s _1 . In fact, 
this mechanism operates most effectively near the dust sub- 
limation radius, 



_ f Lagn 



1/2 



(40) 



1 pc 



/ LAG^ 



V10 46 erg s" 1 



1/2 



1, 200 K 



where <jsb is the Stefan-Boltzmann constant and T su b is 
the dust sublimation temperature. For massive black holes, 
this is within the sphere of influence, where the local escape 
speed is 



v e (R) w 930 km s 



-i / M BH 



VlO 8 M 



1/2 



R 
1 pc 



-1/2 



(41) 



Radiative transfer calculations show that radiation pressure 
on dust in AGN can accelerate winds to comparable ve- 
locities, but not much more (Roth et al. j 2012 1. Since out- 
flows also decelerate substantially as they sweep up ambi- 
ent medium (unless the density profile is very steep), winds 
driven by radiation pressure on dust are therefore limited to 
a few 1,000 km s _1 at most. 

In addition to radiation pressure on electrons and on 
dust, radiation pressure on UV lines (e.g., Murray et al. 



1995 Proga et al.|2000l, hydromagnetic processes (e.g., Em- 


mering et al. 


1992 Konigl & Kartje 19941, and Compton 


heating (e.g., 


Sazonov et al.||2005 1 may contribute to accel- 



crating outflows. An important goal for future work is to 
determine the exact physical mechanisms through which a 
central black hole couples to the surrounding galaxy. 



4.4 Observational signatures of shocked AGN 
wind bubbles 

Energy-conserving outflows predict the existence of shocked 
wind bubbles. We argue here that such bubbles are consis- 
tent with current observations and discuss ways of probing 
them in more detail in the future. 



4-4-1 Direct wind bubble emission 

The free free, inverse Compton, and synchrotron emission 
from the thermal electrons in the shocked wind bubble are 
estimated in Appendix [S] The very hard X-ray (kT c « 
260 keV (T e /3 x 10 9 K)) and relatively faint thermal free 
free will be extremely challenging to detect directly even 
with new-generation observatories such as MiS'LA.Rrl Fur- 
thermore, synchrotron emission from the thermal electrons 
has a characteristic frequency v c ~ 4 Hz j 2 (B/l/iG), which 



http: / /www. nustar.caltech.edu/ 
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is un-observable due to the plasma frequency ~ 3 MHz of 
the Galactic warm interstellar medium | Lacki|2010 1. 

In terms of both luminosity and energy band, the most 
promising signature of the thermal electrons is the inverse 
Compton continuum that they should produce at ~ fcTx ~ 
keV when they interact with the AGN radiation field (see 
also King 2010). At best, when the inverse Compton cooling 
time is short, this signal will be no mor e than a few percent 
of the radiative AGN luminosity (eq. 42 1. Indeed, if Mi n «j n — 
Lagn/c, then a fraction 



(M ia vf n /2) _ Vin 



Lagn 



2c 



= 0.05 



Via 

30, 000 km s" 



(42) 



of the radiative luminosity is converted into wind mechanical 
energy that can then be radiated away. Nevertheless, this 
may be detectable. 



Pounds & Vaughan (2011 1 report a possible detection of 



the inverse Compton continuum in the Seyfert 1 galaxy NGC 
4051, with a luminosity comparable to the mechanical lumi- 
nosity of the wind in that system. Although those authors 
interpret this as evidence of shocked wind cooling and thus 
of a momentum-conserving outflow, this does not necessarily 
follow. In fact, the inverse Compton emission quantifies the 
energy drain from the shocked electrons, but it is the pro- 
tons that provide the majority of the thermal pressure of the 
shocked wind ({j2j. If the electrons experience a significant 
amount of collisionless heating at the shock, but are predom- 
inantly coupled to the protons via Coulomb collisions over 
the longer time scale relevant for shocked wind cooling, then 
inverse Compton emission with a luminosity similar to the 
wind mechanical luminosity could result without removing 
the thermal pressure provided by the protons. 

Relativistic particles accelerated either at the forward 
or reverse shock (e.g., Blandford fc Ostriker|1978 \ may also 
produce non-thermal radio and 7— ray emission. While such 
emission does not provide direct information on the state of 
the gas inside the shocked wind bubble, it probes the me- 
chanical energy processed through the shocks and can con- 
tribute significantly to the spectral energy distribution of 
systems hosting AGN. Kiloparsec-scale radio synchrotron is 



commonly observed from AGN hosts (e.g., Baum et al. 1993 
|Gallimore et al.|2006[ ) and usually interpreted as being pow- 
ered by either a starburst or by a black hole jet interacting 
with the ISM. Cosmic ray electrons accelerated by a wide- 
angle AGN wind may provide an alternative explanation in 
some cases. In particular, Mrk 231 has spatially extended 
radio emission on scales of ~ 100 pc to ~ 25 kpc (Ulves- 
|tad et aLjl 999 ) which the AGN wind can easily account for 
energetically (eq. (B5l). 



4-4-2 Effects of wind bubbles on surroundings 

If the shocked protons do not cool and the outflow is energy- 
conserving, then the hot gas will tend to expand along paths 
of least resistance ({3.21. For an AGN outflow launched at 
the center of a galactic disk, the simplest outcome would be 
bipolar bubbles expanding above and below the disk plane. 



Greene et al. ( 2012 1 report the discovery of such bubbles 



around the luminous obscured quasar SDSS J1356+1026 at 
z — 0.123. The bubbles extend ~ 10 kpc on each side of 
the galaxy, with a probable de-projected expansion veloc- 
ity ~ 1,000 km s _1 . These bubbles, most likely inflated by 



the quasar, are strongly suggestive of buoyant hot gas, such 
as a shocked wind that has not cooled. In this case, the 
observed [OIII] and H/3 emission is consistent with being 
powered by photoionization by an AGN spectrum and thus 
likely traces the swept-up up ambient medium. We must 
note, though, that SDSS J1356+1026 is the only one out of 



15 obscured quasar studied by Greene et al. (2011 1 showing 



such prominent bubbles. Although these bubbles are there- 
fore not ubiquitous in an instantaneous sense, simulations 
of quasars triggered in galaxy mergers indicate that the ac- 
creting black hole may spend up to ~ 90% of its lifetime 
in a buried phase, during which the hot gas may not leak 
effectively. It is thus possible that most of the objects ob- 
served by |Greene et al. ( 2011 1 will eventually inflate bubbles 
similar to SDSS J1356+1026. 

ferm^ also recently detected bubbles of a similar scale 
in our Galaxy, clearly originating from the Galactic Cen- 



ter (Dobler et al. 2010 Su et al. 20101. In our notation, 
the bubbles are inferred to have an age £fl ow = R s /v B ~ 
10 7 yr (Rs/W kpc)(u s /l, 000 km s -1 ) -1 . Interestingly, ob- 
servations also reveal that a starburst in the central par- 
sec occurred 6 ± 2 Myr ago ( Paumard et al.|2006 l, perhaps 
accompanied by an accretion event onto Sgr A* (Levin & 



|Beloborodov 20031. Suppose that the black hole accretion 
event lasted < 1 Myr and inflated a shocked wind bubble 
with T sw ~ 10 9 K (yin ~ 10,000 km s^ 1 ). The bubble may 
have since expanded by a factor ~ 10 in linear dimension, or 
~ 10 3 in volume. Adiabatic cooling (T oc /5 2//3 ) would then 
imply that the bubble would now have an interior temper- 
ature ~ 10 7 K. Thermal free free emission from gas at this 
temperature is consistent with 1.5 keV soft X-ray features 
coincident with the Fermi bubbles previousl y detected by 
ROSAT^\ ( |Bland-Hawthorn fc Cohen] |2003[ ). The bipolar 
bubble morphology would then naturally result from con- 
finement by the massive molecular disk in the inner Galaxy 
(e.g., Zubovas fc Nayakshin|20"T2 1. Thus, the Fermi bubbles 



may be the relics of an energy-conserving AGN outflow. 



4.5 Conclusion and future directions 

Many questions remain open regarding the global impact of 
AGN winds on galaxy evolution. In particular, it is not clear 
if (and if so, how) such winds have a major effect on star 
formation in galaxies. While a quasar outburst can easily 
release enough energy to completely unbind the interstellar 
medium of a galaxy (e.g., Silk fc Rees|1 998), it is not trivial 
to couple this energy effectively to a galactic disk. 

Indeed, galactic disks subtend only a small fraction of 
the solid angle around the black hole and nuclear outflows 
are likely to preferentially escape normal to the disk. Al- 
though OH observations in local ULIRGs suggest short gas 



depletion time scales (Sturm et al. 20111, most of the CO 



flux in some of the same systems is observed to originate 
in a narrow component tracing the star-forming molecular 
disk (e.g., Feruglio et al.|[2010 l, rather than a fast outflow. 



It may be possible to reconcile these observations if the gas 
from the disk is rapidly funneled into the nucleus by strong 
gravitational torques and only gradually incorporated into 



9 http://fermi.gsfc.nasa.gov/ 

10 http:/ /heasarc. gsfc.nasa.gov/docs/rosat/rosat. html 
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the wind. Alternatively, the wind may be gradually mass 
loaded by efficient mixing at the wind-disk interface, ft is 
also possible that the large mass outflow rates have short 
duty cycles, so that they do not necessarily remove a large 
mass fraction. As mentioned above, Greene et al. (20111 
observe significant, non-gravitational disturbances over the 
entire host galaxies of luminous obscured quasars at z < 0.5 
and also find evidence that their interstellar media are com- 
pletely ionized by the quasar. On the other hand, the case 
for effective gas ejection and star formation truncation in 
these systems is less clear. 

The galaxy-scale impact of AGN winds depends on 
three-dimensional effects, non-linear interactions between 



rc Trn = 100 cm 



stellar and black hole feedback processes (e.g., Booth & 
Schaye|2012 1, and the detailed structure of the ISM. It will 
thus be fruitful to incorporate the physical insights gained 
in this work into more detailed numerical simulations. In 
Appendix [C] we outline how our results can be used to im- 
plement the effects of fast nuclear winds more realistically 
in such simulations. 
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APPENDIX A: ENERGY AND MOMENTUM 
CONSERVATION CONDITIONS FOR 
SELF-SIMILAR OUTFLOWS 

In Sj3] we presented examples of numerical solutions relevant 
for observed galaxy-scale AGN outflows and concluded that 
they were energy-conserving. Here we present an analytic 
analysis of self-similar solutions to determine the depen- 
dences on problem parameters and show that energy con- 
servation is realized over a wide parameter space. 

To ensure self-similarity, we assume a power-law ambi- 



ent density profile (eq. (27l) and neglect gravity. We restrict 
ourselves to decelerating solutions with finite free expansion 
radius (see below), corresponding to a < 2. We focus on the 
two-temperature plasma cooling case, which should be re- 
alized for fast nuclear winds in spherical symmetry, and on 
solutions that begin in the energy or momentum-conserving 
limits. We then quantify the transition to a partially radia- 
tive bubble in each case. At this point, the solution devi- 
ates from self-similarity and becomes neither fully energy or 
momentum-conserving. 



Al Solutions that begin energy-conserving 

Consider first solutions that begin in the energy-conserving 
limit. Again assuming that 1/2 of the mechanical energy 
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Figure Al. Radius at which self-similar energy-conserving solu- 
tions transition to a partially radiative bubble as a function of 
the slope of the ambient density profile (solid curves). Fiducial 
AGN parameters (Lagn = 10 46 erg s —1 , r\ n = 1) are assumed, 
and v m = 10,000, 30,000, 50,000 km s" 1 is varied. The dif- 
ferent panels correspond to different densities at the reference 
radius Rq = 100 pc. The free expansion radius, Rf, is shown by a 
dashed curve in each case. The vertical dotted line at a = —1/7 
separates cases where the solution can be energy conserving ev- 
erywhere (when flp RB < Rf for — 1/7 < a < 2) and the cases 
where solutions transition to a partially radiative bubble at a 
finite radius (a < —1/7). 
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Figure A2. Maximum density of the ambient medium at Ro = 
100 pc for which the outflow solution is purely energy-conserving 
(i?p RB < Rf) versus the slope of the ambient medium profile. 
Purely energy-conserving outflows are only possible for — 1/7 < 
a < 2. Fiducial AGN parameters (Lagn = 10 46 erg s —1 , ri n = 1) 
are assumed. 



injected by the AGN goes into kinetic motion of the swept- 
up gas, energy conservation implies 
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L in t = M s v; 
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For a self-similar solutio n of the form R s = Ast^ E ( 



n„„ = l cm 



Rs = PsRs/t), equation ( Al I implies 
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We can then explicitly relate v s and i?. 
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The free expansion radius is defined as the solution to 
M in t = M s (R s (t)): 
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Inside Ri , the integrated mass from the nuclear wind exceeds 
that of the swept-up ambient medium and cooling losses are 
dynamically unimportant. 

In the energy-conserving limit, 
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The transition to a partially radiative bubble occurs when 
this ratio is unity, i.e. at 
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Figure 
eters. 

If and only if —1/7 < a < 2, t p , CO oi/iflow increases with 
R s and solutions can be energy-conserving at all radii past 
R{. Here, the condition a < 2 simply re-iterates the assump- 
tion made earlier to ensure that the wind has a small free 
expansion radius. For fixed AGN parameters and ambient 
profile slope, we define the critical density po,c such that 
Rp RB /Rf = 1. This is the maximum ambient gas density 
at Ro such that the solution is effectively energy-conserving 
everywhere: 
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AGN ' in 



R'n 



Figure 



A2 



shows iij'J 



p§ rit /m p 



as a function of a for 



Ro = 100 pc, fiducial AGN parameters (Lagn = 10 erg 
s _1 , Ti n = 1), and v ia = 10,000, 30,000 and 50,000 km s -1 . 
For example, for a wind with Hi n = 30,000 km s _1 blown 
into an ambient medium of uniform density, the outflow is 
energy-conserving at all times for nn < 8 x 10 6 cm -3 . 

For a < —1/7, solutions can begin energy-conserving 
at Rf but transition to a partially radiative bubble at 



10' 
10 1 

: 10° 

no 1 

10 2 



n H n=100 cm 1 



. 10.000 km/s 

30,000 km/s 

50,000 km/s 
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Figure A3. Radius at which self-similar momentum-conserving 
solutions transition to a partially radiative bubble as a function 
of the slope of the ambient density profile (solid curves). Fiducial 
AGN parameters (Lagn = 10 46 erg s — 1 , r; n = 1) are assumed. 
The different panels correspond to different densities at the refer- 
ence radius Ro = 100 pc. The free expansion radius, R{, is shown 
by a dashed curve in each case. 



Rprb > Rf - However, the conditions for such a transition to 
occur are quite extreme; the large ambient densities required 
would imply much more severe deceleration of the swept-up 
ambient medium than is observed in the systems of interest 
(see the momentum constraint in eq. (|33[)). 



A2 Solutions that begin momentum-conserving 

In a manner analogous to the above, we begin by assuming 
a self-similar momentum-conserving solution R a = Apt 131 " 



(=> v a = /3pRs/t). The equation of motion is then 
Lagn 



yielding 



A P = 



-t = M B Vs 



(3 - a) (4 - a)r in L AG N 



SwpoRqC 



/3p 



and 



v s = 2 



(3 — a)LAGNTin 

8tt(4 - o)cp R% 



1/2 



K 



(a-2)/2 



(P 



(A10) 



(AH) 



(A12) 



(A13) 



While energy-conserving solutions transition to a par- 
tially radiative bubble when ip,cooi/£flow < 1, momentum- 
conserving solutions enter the partially radiative bubble 
phase when t p , CO oi/£cr > 1, where t cr = R sv /vi n . In the 
momentum-conserving phase, _R BW w Rs, so that 



t , 3 8/5 

kp.cool *J 

tcr 8v 2 



a _ 4 \ 3 7v 11 c 19 m 4 e m 1 p i p 6 R 3 a p 3 ov^ 
a-s) e"a«ClnA)*Lj? GN < 



1/10 



(A14) 
x Ri 16 ~ 3a ^ w . 
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This ratio increases with increasing i? s for all a < 2, so 
there is always a transition to a partially radiative bubble 
at a finite radius 



■Rprb 



Figure 



2 35 



-4\ 3 c 19 mXV^ a pg< 
3 J e 16 o"^, (In A) 4 i^GN T in 



l/(3a-16) 



(A15) 



A3 



compares -R PRB with R{ for representative con- 
ditions. The outflow experiences a momentum-conserving 
phase only if -R PRB > Rf, this again requires a fairly high 
ambient density, generally in excess of the mean density 
swept-up by observed galaxy-scale AGN outflows. 



APPENDIX B: EMISSION FROM THE 
SHOCKED WIND BUBBLE 

In this section, we estimate the observational signatures of 
the shocked wind bubble in emission. 

We begin with emission from thermal electrons in the 
shocked wind bubble. For these estimates, we assume the 
late-time, energy-conserving limit. The shocked wind den- 



sity is then given by equation ( 14 1 and the electron temper- 



ature by equation ( 16 1. We neglect the effects of mixing and 
further approximate the bubble volume as 14 = 47r_Rg/3. 
Then, the free free luminosity 

L« ~ esV s 

-1 I 



(Bl) 



4.6 x l(T a erg s" 



Lagn 



\ 10 46 erg s- 



('ii, 



30,000 km s- 



1,000 km s" 



-9/5 



-11/5 



1 kpc 

lnA^ 
40 



11/5 



260 keV (T e /3 x 



with the spectrum peaking around kT c 
10 9 K). 

Similarly, we evaluate the (non-relativistic) inverse 
Compton luminosity, 



Lie ~ 47r 



dRR I j- I C(j T n e (7ph 

y TTIqC J 



(B2) 



9.1 x 10 41 erg s" 1 



La 



10 46 erg s" 



30, 000 km s" 



1,000 km s- 1 



-3/5 



Ra 



-7/5 



1 kpc, 

lnA^ 2/5 
40 



.7/5 



The inverse Compton component should manifest itself at 
the much lower energies characteristic of the AGN radiation 
field, ~ kT x ~ keV ( jKlngpHO L 

The synchrotron emission is sensitive to the unknown 
magnetic field distribution and related to the inverse Comp- 



ton power via equation (19 1. Since the thermal electrons do 



not have ultra-relativistic velocities, the characteristic syn- 
chrotron frequency 



3 7 2 eB 
4n m e c 



4 Hz -y 2 



B 



(B3) 



is so low that the interstellar medium of the Galaxy is 
opaque to this emission | |Lacki|2010"l ). 

Particles can also be accelerated to relativistic veloci- 
ties at the forward and reverse shocks. These can then pro- 
duce non-thermal emission including radio synchrotron, and 
7-rays from inverse Compton scattering and hadronic pro- 
cesses. Although the details are uncertain, a general estimate 
of the non-thermal emission from the AGN shocked wind 
bubble can be obtained by viewing it as a supernova remnant 
analog and exploiting our knowledge of non-thermal emis- 
sion from star-forming galaxies. The gigahertz radio contin- 
uum of star-forming galaxies is in fact believed to originate 
from cosmic ray electrons, and 7— ray emission from cos- 
mic ray protons has also recently been detected in nearby 
starbursts (La cki et al.|2011 |. 

As before, let Li n = Mi u vf n /2 be the mechanical lumi- 
nosity of the AGN wind and denote by E^ NR = /sn/sn,1c£* 
the analogous quantity for SNRs, integrated over the galaxy. 
Here, L* is the stellar bolometric luminosity, /sn is the frac- 
tion going into supernova explosions, and /sN,k is the frac- 
tion of the supernova energy going into kinetic power (rather 
than being radiated away). Then: 



L in 



5Ti 



Lagn 
L* 



0.01 



(B4) 



/sN,k 



Lagn 



where the fiducial value for /sn is based on the results of 



Abbott (19821. Assuming a comparable efficiency for accel- 



erating particles in SNR and AGN wind shocks, this shows 
that the non-thermal emission from the AGN wind bubble 
can in principle significantly exceed that powered by star 
formation. In a ULIRG like Mrk 231 with a powerful quasar 
at the center, Lagn ~ L+ (e.g., Veilleux et al. 2009). We cau- 
tion, however, that this ratio should be taken with a grain 
of salt since the conditions near the shock acceleration sites 
in SNRs and in the shocked wind bubble (including the pho- 
ton and magnetic energy densities, and the time scales for 
advection and adiabatic losses) need not be the same. 

Of particular interest is the non-thermal synchrotron 
contribution, which may have been observed in some sys- 
tems. An upper bound can derived under the assumption 
that the system acts as a calorimeter, i.e. that the cosmic 
ray electrons radiate all their energy in synchrotron on a 
time scale short relative to all other radiative and dynamical 
time scales. Assuming further that the particles are shock 
accelerated with a spectrum 71(7) oc 7 _p with p = 2 (Bland- 



ford & Eichler 19871, the synchrotron emission has equal 



power per frequency decade, 



Lmax 
v\ \. 
" I synch 



iLm _ ^T in 

ln7 ma x 2hi7 n 
3.6 x 10 41 erg s _1 t- u 



(v) Lagn 



(B5) 



£ \ /In 7 



0.01 ) V InlO 6 

fin \ / -I/AGN 



0.1c/ V10 46 erg s- 1 



The efficiency factor £ ~ 1% is constrained for SNR shocks 
based on modeling the far infrared-radio correlation (e.g., 
Lacki et al.|2010 |. The synchrotron luminosity estimated in 
equation (B5l can only be realized if Ub > U p h, which may 
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not be realized given the large photon densities near AGN 



(see eq. (18l) 



APPENDIX C: MODELING THE EFFECTS OF 
FAST NUCLEAR WINDS IN SIMULATIONS 

It is at present extremely challenging to simultaneously re- 
solve the launching of AGN winds and their effects on galaxy 
scales in numerical simulations. For example, the accretion 



disk wind model of Murray et al. ( 1995 ) predicts that winds 
are launched at ~ 0.01 pc from the black hole. However, one 
can model their effects in galaxy simulations by imparting 
kinetic kicks to gas elements, and attempt to resolve the 
resulting wind bubble (e.g. 



DeBuhr et al. 2011a I. In this 



section, we use our analytic model to outline the basic re- 
quirements for numerical simulations to faithfully capture 
the impact of fast nuclear winds on galaxy scales. It is not 
our intention to provide a complete algorithmic description, 
as this would depend on the details of the numerical scheme 
used, but rather to summarize the criteria that realistic im- 
plementations must achieve. 

As discussed in it is the cooling of the shocked wind, 
rather than the shocked ambient medium, that determines 
whether the outflow conserves energy or momentum. Thus, 
a realistic simulation should at minimum resolve the wind 
shock. The mass of the shocked wind is 

Rs 



Ms- 



2.2 x 10 6 Mq / e 



(CI) 



1 kpc 



0.1 

000 km s- 



V 10 8 Mq 

-1 



In the energy-conserving phase, the radius of the wind shock 
is 

1/2 



Rsw ^ Rs 



t'h- 



(C2) 



180 pc 



1 kpc J \1, 000 km s" 1 

1/2 



1/2 



30,000 km s- 1 



whereas i? S w ~ Rs when it is momentum-conserving. How- 
ever, a significantly better resolution than suggested by the 
above estimates is needed in general to accurately capture 
the cooling time of the shocked wind. In particular, many 
numerical codes broaden shocks over several resolution ele- 
ments, which allows gas to artificially cool while it is being 



shocked (Hutchings & Thomas 2000 Creasey et al. 20111 



As a result, the simulated gas may never reach the physical 
post-shock temperature 



10 10 K 



Uin/30,000 km s -1 ) 5 



(C3) 



Thus, numerical simulations should: 



• Ensure that a wind shock is present, with radius, 
shocked wind mass and temperature consistent with equa- 
tions ( C1|C3 1 in the energy-conserving phase. 

• Implement a cooling function extending to T > 10 10 
K, including Compton scattering. Two-temperature plasma 



effects in the shocked wind can be modeled following the 
results of ij2] The cooling time of the shocked wind should 
be checked against the analytic expectations to diagnose po- 
tential artificial losses due to numerical effects. 

• The time integration scheme must ensure accurate en- 
ergy conservation at the numerical level. For instance, in 
the absence of a proper limiter neighboring smooth parti- 
cle hydrodynamics particles that are initially at rest may 
not react adequately to sudden changes, which can lead to 
effects including inter-particle crossing and violation of en- 



ergy conservation (Saitoh & Makino 2009 Durier & Dalla 
Vecchia|2012 l. 



If it is not possible to correctly resolve the dynamics of the 
wind bubble, approximations that enforce physical energy 
conservation may be considered. For instance, some imple- 
mentations of supernova feedback suppress cooling in order 
to avoid efficient energy radiation when the ISM is not prop- 
erly resolved (for a review, see Kay et al. 2002; Dalla Vecchia 
& Schaye 2012 discuss alternatives), and other techniques 
have also been developed in the context of AGN feedback 
(e.g., |Sijac"ki et al.|200f| |Booth fc Schaye|2009| >. 
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